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Abstract. In this paper we present and analyse determinations of effective temperatures of planet-hosting stars using infrared 
(IR) photometry. One of our goals is the comparison with spectroscopic temperatures to evaluate the presence of systematic 
effects that could alter the determination of metal abundances. To estimate the stellar temperatures we have followed a new 
approach based on fitting the observed 2MASS IR photometry with accurately calibrated synthetic photometry. Special care has 
been put in evaluating all sources of possible errors and incorporating them in the analysis. A comparison of our temperature 
determinations with spectroscopic temperatures published by different groups reveals the presence of no systematic trends 
and a scatter compatible with the quoted uncertainties of 0.5-1.3%. This mutual agreement strengthens the results of both the 
spectroscopic and IR photometry analyses. Comparisons with other photometric temperature calibrations, generally with poorer 
performances, are also presented. In addition, the method employed of fitting IR photometry naturally yields determinations 
of the stellar semi-angular diameters, which, when combined with the distances, results in estimations of the stellar radii with 
remarkable accuracies of ~2— 4%. A comparison with the only star in the sample with an empirically determined radius (HD 
209458 - from transit photometry) indicates excellent agreement. 
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1. Introduction 

The characterization of the properties of planet-hosting stars 
has been an active field of study. Soon after the discovery of the 
first candidates, claims were made that stars with planets dis- 
played on average higher metal contents (Gonzalez 1997) than 
other solar neighbourhood stars. A number of subsequent in- 
dependent studies with increasingly large stellar samples have 
mostly confirmed the initial claims (e.g. Santos et al. 2003 ). An 
important point to be made is that the determination of chem- 
ical abundances, mostly carried out through detailed analysis 
of spectroscopic data, is quite challenging (see Gonzalez 2003 
for a complete review). As it has been shown for late-type stars 
(i.e. FGK planet hosts), a strong degeneracy affecting the deter- 
mination of metal abundances is the correlation with effective 
temperature. Straightforward estimations show that a system- 
atic error of +100 K in T e ff (i.e. 1.5-2% at the temperatures 
of FGK stars) results in metal abundances being systematically 
overestimated by +0.06 dex (~15%). Most studies of stellar 
atmospheric parameters carry out multiple fits to derive chem- 
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ical compositions and effective temperatures from the spectra. 
Although the aforementioned correlation would not alter the 
conclusions of relative studies including both planet-hosting 
and non-planet hosting stars, the metal richness of the stars 
would be systematically biased when compared to the Sun. 
Another point worth making is the use by present spectroscopic 
studies of solar line oscillator strengths for all late-type stars 
(e.g. Gonzalez & Laws 2000 1. This might introduce systematic 
errors for temperatures below and above that of the Sun that 
thus far have not been addressed in detail. 

The potential problems with spectroscopic analyses dis- 
cussed above make a completely independent temperature de- 
termination, for example using photometry, very valuable. The 
absence of systematic effects when comparing photometric and 
spectroscopic temperatures would strengthen the case for the 
metal richness of planet-bearing stars and support the use of 
solar oscillator strengths over the relevant spectral type range. 
However, the determination of photometric temperatures for 
cool stars (below 7000 K) is not straightforward because most 
photometric systems are not designed for such low tempera- 
tures. For example, although some efforts have been made to 
extend the temperature range covered by Stromgren calibra- 
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tions down to late-type stars (Olsen 1984 1, most of the work is 
still in a preliminary stage. Here we present a new approach, 
namely the determination of effective temperatures from in- 
frared (IR) photometry. The underlying idea is similar to the 
Infra-Red Flux Method (IRFM), proposed and implemented by 
Blackwell & Shallis fT9T7l . Blackwell et al. (1980 1 and later 
Alonso et al. (I1996> . 

In this paper we briefly discuss the proposed approach and 
compare our results with spectroscopic temperature determina- 
tions. In addition, the analysis also yields an robust and accu- 
rate determination of the stellar radius, provided the distance 
is known. With the release of the 2MASS All Sky Catalog 1 , 
which contains IR photometry covering the entire sky, the pro- 
posed method has a wide applicability, thus permitting accurate 
(a few percent) and effortless determinations of temperatures 
and radii (important for transits) of planet-hosting stars. 

2. The stellar sample 

Planet searching projects are currently so efficient that the num- 
ber of planet-hosting stars increases on a weekly basis. Our 
starting sample was complete as of July 2003 and comprised 
94 stars with planets, which makes up a statistically signifi- 
cant number for our comparisons. The original list was cross- 
matched with the 2MASS All Sky Catalog with a 100% success 
rate. In a next step, we excluded those stars with poor-quality 
photometry in all of the three 2MASS bands, generally because 
of strong saturation. 2 As a criterion, we rejected stars with er- 
rors in the K band (the band with best quality for bright stars) 
greater than 0.05 mag. In addition, two stars were excluded for 
different reasons: HD 113020, classified as spectral type M4, 
because its temperature is expected to be below 4000 K and 
beyond the studied temperature range; and OGLE-TR56, the 
farthest known planet-hosting star at 1 .5 kpc, because of con- 
cerns with the extinction correction. The final working sam- 
ple is therefore composed of 8 1 planet-bearing stars, which are 
listed in Tabled 

The 2MASS photometry was complemented with V stan- 
dard measurements (together with their uncertainties) obtained 
from GCPD (Mermilliod et al. [T997t . Also from GCPD we 
compiled Stromgren photometry and carried out an estimation 
of the reddening correction. All stars were found to fall right 
on the standard relationships as defined by Crawford i 19751 1 
and Olsen ( 119841 IT988 ). i.e. zero reddening as expected from 
their nearby distances within 100 pc. 

3. Determination of temperatures and radii 

The use of IR photometry for determining effective tempera- 
tures was initially proposed by Blackwell & Shallis ( 1977 1. 
Their so-called IRFM uses the ratio between the bolometric 
flux of the star and the monochromatic at a given wavelength 

1 http://www.ipac.caltech.edu/2mass 

2 The 2MASS team managed to extract good quality photometry 
from the wings of mildly saturated stars, which led to a significant 
increase in dynamic range (see Cutri et al. 2003 l. In general, good 
photometry was obtained for stars with K > 4 mag. 



in the IR, both measured at the Earth, as the observable quan- 
tity. This ratio is compared with a theoretical estimate derived 
from atmosphere models, allowing the determination of the ef- 
fective temperature. The method we have developed follows a 
somewhat different approach, namely, a fit of the stellar spec- 
tral energy distribution from the optical (V) to the IR (JHK) 
with synthetic photometry computed from atmosphere models. 
This wavelength range ensures enough sensitivity to temper- 
ature, while being essentially immune to intrinsic variations 
(mostly due to chromospheric activity) that would strongly af- 
fect bluer passbands. Also importantly, the spectral energy dis- 
tribution in the optical/IR for cool stars is mostly insensitive to 
the (often unknown) values of log g and [Fe/H]. In contrast to 
the IRFM, our method does not require the use of a bolometric 
flux calibration, thus reducing the risk of uncontrolled system- 
atic errors, and permits the determination of individual uncer- 
tainties. The algorithm will be discussed in detail by Masana 
(2004 1 but we briefly outline the procedure here. 

A critical pre-requisite to compare stellar and model fluxes 
is a well-characterized photometric system and an accurate 
flux calibration. Fortunately, the recent work by Cohen et al. 
(2003a 2003b I provides a set of consistent absolute flux cali- 
brations in both the optical (Landolt system) and IR (2MASS). 
To compute synthetic magnitudes we employ the dense grid of 
Kurucz atmosphere models 3 , although tests with other models 
(discussed below) reveal negligible differences. The synthetic 
VJHK photometry was calculated for each of the grid points as 
a function of T e g, [Fe/H], and logg. We restricted our calcu- 
lations to the temperature interval between 4000 and 8000 K. 
The upper limit is defined by the increasing dependence of the 
results on log g and the lower limit is set by the decreasing per- 
formance of the models because of molecular bands. 

The algorithm itself is based on the minimisation of a 
X 1 function defined from the difference between the observed 
and synthetic VJHK magnitudes. The two adjustable parame- 
ters are T e s and a magnitude zeropoint (ZP), while \ogg and 
[Fe/H] are fixed input data (recall that interstellar absorption 
for our stars is negligible). The ZP is the difference between 
the synthetic (star's surface) and the observed magnitude (Earth 
surface), and it is directly related to the semi-angular diameter 
(0 = 10" a2zp ). The input values of [Fe/H] were taken from 
the spectroscopic determinations and logg was roughly esti- 
mated from stellar models (in the T^-L plane) after one itera- 
tion. Convergence towards the minimum value of the reduced 
X 2 was very fast and reached values close to unity. 

Perhaps the greatest advantage of our method is that it 
yields the individual uncertainties of both T s $ and 8, calcu- 
lated from the covariance matrix, provided the involved errors 
are realistically estimated. For the present sample, the error 
sources considered are: 1) Observed magnitude errors, with an 
minimum 0.015 mag error in V (usually computed from the 
average of a few measurement and thus prone to error under- 
estimation); 2) Zero point uncertainties in the flux calibration 
(1.5-1.7%) from Cohen et al. ( I2003all2003bl : 3) Conservative 



3 http://kurucz.harvard.edu/grids.html 
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error bars of 0.1 dex in [Fe/H] and 0.2 dex in logg dex 4 ; 4) 
No error was attributed to the model fluxes. We ran a number 
of comparison with other (less dense) atmosphere model grids 
such as those by Castelli et al. ( 1997 1 and the NextGen models 
by Hauschildt et al. 119991 . Our tests were extremely satisfac- 
tory and yielded average temperature differences below 15 K 
(~0.3%) in all cases. Thus, systematic errors introduced by at- 
mosphere models are likely to be negligible. 

The procedure described above yields two basic parame- 
ters: The best-fitting effective temperature and semi-angular di- 
ameter (i.e. radius/distance). With a known distance, the latter 
can be transformed into a true radius measurement. The results 
for our stellar sample are listed in Tabled The relative accu- 
racy of the results is 0.5-1 .3% in effective temperature and 0.9- 
2.4% in semi-angular diameter. Since the planet-hosting stars 
in the sample are generally nearby, their Hipparcos distances 
are very accurate. This results in relative uncertainties in the 
individual radii in the range of \3-A.l% for 90% of the stars. 

4. Discussion 

With the temperatures from IR photometry in hand, we carried 
out a comparison with spectroscopic determinations to assess 
their mutual agreement. We considered two comprehensive and 
independent spectroscopic analyses, following, however, simi- 
lar approaches: The Swiss group (Santos et al. 2001 , 2003 1 and 
the U S group (Gonza lez [T9971 1 1 9981 1 1 9991 Gonz alez & Laws 
120001 Gonzalez et al. l200Tl Feltzing & Gonzalez I200T1 Laws 
et al. l2003> . A cross-match of the IR and spectroscopic sam- 
ples yields a total of 69 stars in common with the Swiss group 
and 49 with the US group. The resulting comparisons are illus- 
trated in Fig.^ where the observational data are plotted with 
their error bars. As can be seen, the agreement among the three 
temperature determinations is very remarkable. The mean dif- 
ference between the IR-based and the Swiss group tempera- 
tures is < r e ff IR - r e ff sw >= _ 4.2 + 5 1 .4 K, with no hint of any 
systematic trends neither as a function of temperature, metallic - 
ity or surface gravity. The difference with the US group yields 
a value of < r e ff IR - T e ff US >= +20.5 + 65.9 K, with a low 
significance trend of larger IR temperatures at the low end 
and smaller IR temperatures at the high end. No systematic 
trends are seen as a function of metallicity or surface gravity. 
Interestingly, in both cases the scatter of the differences is en- 
tirely consistent with the quoted error bars, which indicates that 
our procedure yields realistic uncertainties. 

A valuable independent test of our analysis can be carried 
out by comparing the semi-angular diameters with observa- 
tional determinations. Unfortunately, this is not possible at this 
point because all the stars with available empirical angular di- 
ameter measurements are very bright (see Alonso et al. 1994 
for a compilation) and do not have accurate 2MASS photom- 
etry. Luckily, there is one star in our sample with an accurate 
empirical radius measurement and this is HD 209458, which 
undergoes planetary transits. Generally it is very difficult to 

4 Our method is essentially insensitive to the adopted value of 
[Fe/H] and logg. Thus, uncertainties up to 0.3 dex in [Fe/H] and 
0.5 dex in logg introduce uncertainties in r eff below 0.5%. 
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Fig. 1. Comparison of effective temperatures determined us- 
ing our method (IR) and temperatures determined spectroscop- 
ically by the Swiss (SW) and US groups (see text). 

disentangle the stellar and the planetary radii effects when ana- 
lyzing transit light curves. However, Brown et al. (2001 1 were 
able to determine both the stellar radius and the planetary ra- 
dius simultaneously from a high precision transit light curve of 
HD 209458 obtained with the Hubble Space Telescope. Their 
empirically-determined radius is 1.146 + 0.050 R©, which can 
be compared with our estimate of 1.145 + 0.049 R from the 
fit to the star's optical/IR energy distribution. The agreement 
is excellent (owing the almost null difference to small number 
statistics) and the error bars of both measurements are alike. 

These important results indicate that our temperature and 
radius determination method is robust and with a very similar 
accuracy to the one possibly achievable with the best available 
spectroscopy and photometry today. 

For completeness, the effective temperatures in Table ^ 
were also compared with those computed from photometric 
calibrations. A detailed discussion is left for a forthcoming 
publication, but we shall briefly review here some of the most 
relevant results. Our tests focused on the Stromgren calibra- 
tions of Olsen ( 1984 1 and the IRFM as implemented by Alonso 
et al. (I1996> . The comparisons show that the (b - y) calibra- 
tion of Olsen ( 1984 1 yields effective temperatures that are sys- 
tematically lower than those resulting from our method, with 
the mean difference being 168+103 K. Interestingly, a marked 
trend was found when plotting the differences as a function of 
metallicity. It seems that the origin of the discrepancies can be 
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Table 1. Effective temperatures, semi-angular diameters (6 = radius/distance) and radii, with their uncertainties, of planet-hosting 
stars as resulting from the algorithm discussed in this paper. 
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ascribed to the stellar sample used to derive the Stromgren cal- 
ibrations, since only one third of Olsen's reference stars have 
super-solar metallicities. In contrast, the comparison with the 
temperatures derived from the IRFM calibration of Alonso et 
al. {1996) reveals no systematic trends but they are on average 
76+42 K smaller than our determinations. Given the similarity 
of both methods and input data, the discrepancy can tentatively 
be attributed to the photometric transformations employed or 
to the bolometric flux calibration used by Alonso et al. (I1996> . 

5. Conclusions 

The conclusions of our study are twofold. First, we have 
compared the effective temperature determinations for planet- 
bearing stars from two completely independent approaches 
with similar accuracies, namely detailed spectroscopic analy- 
sis and IR photometry. The results indicate an excellent agree- 
ment in the entire temperature range, which confidently rules 
out the possibility of systematic effects in spectroscopic metal- 
licity determinations and supports the use of solar line oscil- 
lator strengths. Second, the method presented, consisting in a 
fit to the observed VJHK magnitudes using synthetic magni- 
tudes, has proved its reliability, yielding accurate (~1%) and 
cost-effective temperatures. As a bonus, the analysis also pro- 
vides determinations of the semi-angular diameters and, even- 
tually, the stellar radii. The resulting radius accuracy of a few 
percent (for nearby stars) could be extremely useful to break 
the strong degeneracy between the radii of the planet and the 
star when analysing transit light curves. 
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